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1  Overview 

The  connections  between  dissipation  inequalities  and  integral  quadratic  constraints  (IQCs)  was  the  major 
thrust  explored  during  the  program.  In  particular,  it  was  shown  that  existing  frequency  domain  IQC 
analysis  conditions  are  equivalent  (under  mild  technical  assumptions)  to  a  related  time-domain  dissipation 
inequality  condition.  This  time-domain  approach  enables  applications  of  IQCs  to  analyze  the  robustness 
of  uncertain  time-varying  and  nonlinear  systems.  These  theoretical  techniques  can  be  used  to  improve  the 
design  and  robustness  of  advanced  flight  control  algorithms. 

This  report  documents  the  research  performed  as  part  of  the  project  entitled  “Development  of  Nonlinear 
Analysis  Tools  Based  on  a  Merged  IQC/SOS  Theory”.  This  research  is  funded  by  the  AFOSR  under  grant 
FA9550-12-1-0339.  The  technical  monitor  for  the  program  is  Dr.  Fariba  Fahroo.  The  following  subsections 
summarize  the  key  findings  of  the  research  supported  on  this  contract.  Details  can  be  found  in  the 
publications  listed  in  Section  7. 

2  Dissipation  Inequalities  and  Integral  Quadratic  Constraints 

Integral  quadratic  constraints  (IQCs),  introduced  in  [9-11],  provide  a  general  framework  for  robustness 
analysis.  In  this  framework  the  system  is  separated  into  a  feedback  connection  of  a  known  linear  time- 
invariant  (LTI)  system  and  a  perturbation  whose  input-output  behavior  is  described  by  an  IQC.  The 
IQC  stability  theorem  in  [9-11]  was  formulated  with  frequency  domain  conditions  and  was  proved  using 
a  homotopy  method.  The  remainder  of  this  section  briefly  describes  stability  theorems  using  dissipa¬ 
tion  inequalities  and  integral  quadratic  constraints.  The  main  contribution  of  the  work  was  to  show  an 
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equivalence  between  dissipation  theory  and  IQC  approaches  [8,21,23].  The  benefit  of  the  time-domain  dis¬ 
sipation  inequality  approach  is  that  it  can  be  generalized  to  cases  where  the  known  plant  in  the  feedback 
interconnection  is  nonlinear  and/or  time- varying.  For  example,  dissipation  inequality  conditions  for  linear 
parameter  varying  systems  [24]  can  be  extended  to  include  uncertainty.  Details  on  this  work  can  be  found 
in  Reference  2  of  the  publications  listed  in  Section  7. 

2.1  Problem  Formulation 

Consider  the  feedback  interconnection  shown  in  Figure  1.  This  interconnection  is  specified  by  the  following 
equations: 

v  —  Gu  T  fj  u  =  A(v)  +  r  (1) 

where  r  G  L^efO,  °o)  and  /  G  Z/fJC/oo)  are  exogenous  inputs.  A  :  L^JO,  oo)  -Q^Cfoo)  is  a  causal 
operator  with  bounded  gain.  G  is  a  linear  time-invariant  system: 

xq  =  Axq  +  Bu,  y  =  Cxq  +  Du  (2) 

where  xq  G  Kn<3  is  the  state  of  G. 


Figure  1:  Feedback  interconnection 

Definition  1  The  interconnection  of  G  and  A  is  well-posed  if  for  each  r  G  L'fff  0,oo)  and  f  G  L^e  [0,  oo) 
there  exist  unique  u  G  L^[0, oo)  and  v  G  such  that  the  mapping  from  (?’,/)  to  (u,v)  is  causal. 

Definition  2  The  interconnection  of  G  and  A  is  stable  if  it  is  well-posed  and  if  the  mapping  from  (r,  /) 
to  (u,v)  has  finite  L2  gain  for  all  solutions  starting  from  xq( 0)  =  0. 

2.2  Frequency  Domain  IQC  Stability  Condition 

Let  II  :  jWL  — »•  <[An+m) x  (n+m)  ^g  a  measurable  Hermitian- valued  function.  Two  signals  v  G  L^/O,  °o)  and 
w  G  L™[0,  00)  satisfy  the  IQC  defined  by  the  multiplier  II  if 

/oo  r  -1  *  1-  1 

3S  du  >0  (3) 

-00  L  J  L  J 

where  v(juj)  and  w(ju>)  are  Fourier  transforms  of  v  and  w.  A  bounded,  causal  operator  A  :  L^fO,  00)  — > 
Lffl  0, 00)  satisfies  the  IQC  defined  by  II  if  Equation  3  holds  for  all  v  G  [0,  00)  and  w  =  A(u).  The  next 
theorem  provides  a  stability  condition  for  the  interconnection  of  G  and  A. 

Theorem  1  (  [11])  Let  G  G  Ki^™  and  A  :  Ufe  — >  Lff  be  a  bounded  causal  operator.  Assume  for  all 
t  G  [0, 1]: 

1.  the  interconnection  of  G  and  r  A  is  well-posed. 

2.  t A  satisfies  the  IQC  defined  by  II. 
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3.  Be  >  0  such  that 


G(j^  1  *  n (ju)  \  G^u)  1  <  -el  Vw  €  R.  (4) 

Then  the  feedback  interconnection  of  G  and  A  is  stable. 

For  rational  multipliers,  Condition  3  is  equivalent  to  an  LMI.  Specifically,  any  II  G  x  Gl+mi 

can  be  factorized  as  II  =  where  M  =  MT  G  WT'zXnz  and  \F  G  WH^x^n+m\  Such  factorizations 

are  not  unique  but  can  be  computed  with  state-space  methods  [16].  Denote  a  state-space  realization  of 
’P  by  (A^,  [B^i,  B^2],  C^,  [D^i,  D^2])  where  the  B^/D^,  matrices  are  partitioned  compatibly  with  [^].  A 
state-space  realization  for  the  system  \F  [  cj  ]  is: 

(A,  B,  C,  D)  :=  ^  B^c  At  ,  B^+B^D  ,[D-4’iC  Cy,]  ,  D^2  +  A/u-D)  (5) 

Finally,  the  KYP  Lemma  [14,20]  can  be  applied  to  demonstrate  the  equivalence  of  Condition  3  in  Theorem  1 
to  an  LMI  condition.  This  result  is  stated  formally  below. 


Theorem  2  3e  >  0  such  that  Equation  4  holds  if  and  only  if  there  exists  a  matrix  P 


ATP  +  PA  PB\  \CT' 
BtP  0  J  +  [dt_ 


M  [C  D]  <  0 


PT  such  that 


(6) 


2.3  Time  Domain  Dissipation  Inequality  Stability  Condition 

An  alternative  time-domain  stability  condition  can  be  constructed  using  IQCs  and  dissipation  theory.  Let 
(\F,M)  be  a  factorization  of  II.  Let  signals  (v,w)  satisfy  the  IQC  in  Equation  3  and  define  z(ju)  := 
'F(jw)  ■  Then  the  IQC  can  be  written  as  JfQ  z(jcc)* Mz(ju)duj  >  0.  By  Parseval’s  theorem  [27], 

this  frequency-domain  inequality  can  be  equivalently  expressed  in  the  time-domain  as: 

z(t)T  Mz{t)  dt  >  0  (7) 

where  z  is  the  output  of  the  LTI  system  T: 

ip(t)  =  A^ip(t)  +  B^iv(t)  +  B^2w(t),  if(0)  =  0  (8) 

z(t)  =  C^(t)  +  D^iv(t)  +  D^2w(t)  (9) 

Thus  A  satisfies  the  IQC  defined  by  II  =  \P~M\P  if  and  only  if  the  filtered  signal  z  =  T  [  ^  ]  satisfies  the 
time  domain  constraint  (Equation  7)  for  all  v  G  Ufl 0,  oo)  and  w  =  A(u). 

The  constraint  in  Equation  7  holds,  in  general,  only  over  infinite  time.  The  term  hard  IQC  in  [11]  refers 

to  the  more  restrictive  property:  z(t)T Mz(t)  dt  >  0  holds  VT  >  0.  In  contrast,  IQCs  for  which  the 

time  domain  constraint  need  not  hold  for  all  finite  times  are  called  soft  IQCs.  This  distinction  is  important 
because  the  dissipation  theorem  below  requires  the  use  of  hard  IQCs.  One  issue  is  that  the  factorization 
of  II  is  not  unique.  Thus  the  hard/soft  property  is  not  inherent  to  the  multiplier  II  but  instead  depends 
on  the  factorization  (\F,M).  A  more  precise  definition  is  now  given. 

Definition  3  Let  II  G  be  factorized  as  \P~ M \F  where  M  =  MT  G  MnzXTlz  and  T  G 

MHto  x(n+m^ .  Then  (T,  M )  is  a  hard  IQC  factorization  of  II  if  for  any  bounded,  causal  operator  A  satisfying 

the  IQC  defined  by  II  the  following  inequality  holds 

z{t)T M z(t)  dt  >  0  (10) 

for  all  T  >  0,  v  G  L^O, oo),  w  =  A(v),  and  z  =  \F  [£,]. 
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The  stability  of  the  feedback  system  can  be  analyzed  using  Figure  2.  This  feedback  interconnection 
including  T  is  described  by  w  =  A(v)  and  the  following  extended  linear  dynamics  (omitting  the  dependence 
of  all  signals  on  time  t): 

x  =  Ax  +  Bw  +  B2  [(]:=  F(x,w,f,r)  (11) 

[vu]  =  C1x  +  D11w  +  D12[f]  (12) 

z  =  Cx  +  Dw  +  t)22  [{]  (13) 


where  x  :=  [xq,  bT]T  £  is  the  extended  state.  A,  B,  C ,  and  D  are  defined  in  Equation  5.  The 

remaining  state  matrices  are  defined  as: 


B2  := 


0  B 
l  x  D 


,  Ci  :  = 


(14) 


Du  :  =  [  ]  ,  D12  ■=  [o'/])  b22  :=  [  D^D  ] 


(15) 


Figure  2:  Analysis  Interconnection  Structure 

The  next  theorem  provides  a  stability  condition  using  IQCs  and  a  standard  dissipation  argument. 

Theorem  3  Let  G  £  li”xm  and  A  :  LV2e  — y  L™e  be  a  bounded  causal  operator.  Assume  that: 

1.  the  interconnection  of  G  and  A  is  well-posed. 

2.  A  satisfies  the  IQC  defined  by  II  and  (tH ,M )  is  a  hard  factorization  of  II. 

3.  there  exists  P  >  0  and  a  scalar  7  >  0  such  that  V{x)  :=  xTPx  satisfies 

ztMz  +  VI/  •  F(x,w,  f,r)<  7  [/f  [/]  -  ^  [vf  [“]  (16) 

for  all  nontrivial  (x,  w,  r,  f)  £  Mnc+r7’  x  Mm  x  Mm  x  Mn  where  u,  v,  z  are  defined  by  Equations  12  and 
13. 

Then  the  feedback  interconnection  of  G  and  A  is  stable. 


Equation  16  is  an  algebraic  inequality  on  the  variables  ( x,w,f,r ).  This  constraint,  when  evaluated 
along  solutions  of  the  extended  system,  represents  the  differential  form  for  a  dissipation  inequality  satisfied 
by  the  extended  system.  The  next  lemma  shows  that  the  dissipation  inequality  in  Equation  16  is  also 
equivalent  to  the  KYP  LMI. 


Lemma  1  There  exists  P  >  0  satisfying  the  dissipation  inequality  (Equation  16)  for  some  7  >  0  if  and 
only  if  there  exists  P  >  0  satisfying  the  KYP  LMI  (Equation  6). 
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2.4  Main  Result 


The  previous  section  summarizes  two  IQC  stability  theorems.  Theorem  1  involves  a  frequency  domain 
condition  with  a  multiplier  II.  Theorem  3  involves  a  dissipation  inequality  with  a  multiplier  (T,  M).  The 
multipliers  are  connected  by  a  non-unique  factorization  II  =  Theorems  1  and  3  are  clearly  related 

as  the  Condition  3  in  each  theorem  is  equivalent  to  the  same  KYP  LMI.  Two  important  properties  are 
required  for  the  dissipation  inequality  approach: 

1.  (4 f,M )  must  be  a  “hard”  factorization  to  ensure  the  time-domain  constraint  holds  over  all  finite 
intervals. 

2.  The  solution  to  the  KYP  LMI  must  satisfy  P  >  0.  This  is  not  required  for  the  frequency  domain 
test. 

The  main  result  is:  For  a  class  of  multipliers,  II  has  a  factorization  (T,  M)  that  is  both  “hard”  and  such 
that  any  feasible  solution  of  the  KYP  LMI  satisfies  P  >  0. 

2.4.1  Condition  for  Hard  Factorization 

Define  the  following  cost  functional  J  on  v  G  LI];  [0,  oo),  vj  G  L™[ 0,  oo),  and  t/’o  G  MnY 

poo 

J(v,  w,  V’o )  :=  /  z(t)TMz(t)dt  (17) 

Jo 

subject  to: 

i>(t)  =  A^ip(t)  +  B^v{t)  +  B^w(t),  if(0)  =  V’o 
z(t)  =  +  A/, iv(f)  -I-  D^2w(t) 


Also  define  the  upper  value  J  as 

J{ip o)  :=  inf  sup  J(v,w,tf o)  (18) 

veL%[0,oo)  weL?\o,oo) 

Lemma  2  Let  II  G  x  (n+m)  he  a  multiplier  and  (4 r,M)  any  factorization  of  U  with  T  stable. 

Assume  A  :  L^O,  oo)  — >  L2g[0,oo)  is  a  casual,  bounded  operator  that  satisfies  the  IQC  defined  by  II.  Then 
for  all  T  >  0,  v  G  L’2”[0,  oo)  and  w  =  A(u),  the  output  o/T  satisfies: 

z(t)T M z(t)  dt  >  —  J (V’t)  (19) 

where  ifx  denotes  the  state  of  T  at  time  T  when  driven  by  inputs  ( v,w )  with  initial  condition  if(0)  =  0. 

2.4.2  Condition  for  Positive  Semidefinite  KYP  Solution 

Define  the  lower  value  J  as 


Jfip o)  :=  sup  inf  J(v,w,ip o) 

■weL™  [0,oo)^eL^[0,oo) 


(20) 


Lemma  3  Let  II  G  ML^+m'|x('n+m'>  be  a  multiplier  and  (\H,M)  any  factorization  of  II  with  T  stable.  Given 
G  G  assume  the  corresponding  KYP  LMI  (Equation  6)  is  feasible  with  state  matrices  ( A ,  B ,  C,  D) 

defined  in  Equation  5.  Let  P  =  PT  denote  a  solution  to  the  KYP  LMI.  Then  K(.xo)  :=  XqPxq  >  J(V’o) 
for  all  s0  :=  [  *5,0. 1>o  }T  G  RnG+n^ . 
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2.4.3  Dissipation  Inequalities  with  J-Spectral  Factorizations 

By  Lemma  2,  (4',  M )  is  a  hard  factorization  if  <  0  V^.  By  Lemma  3,  all  KYP  LMI  solutions  satisfy 
P  >  0  if  >  0  V'i/’.  Moreover,  weak  duality  implies  that  the  lower  and  upper  values  satisfy 
Hence  a  factorization  n  =  vH~Mvl/  that  is  both  “hard”  and  ensures  P  >  0  for  all  KYP  LMI  solutions  must 
have  0  <  )  <  J(V’)  £  0.  In  other  words,  for  such  a  factorization  the  lower  and  upper  values  must 

satisfy  =  J(ip)  =  0.  The  following  special  factorization  plays  a  key  role  in  the  main  result  below. 

Definition  4  (1L,M)  is  called  a  Jn,m- spectral  factor  of  n  =  n~  £  ML^+m^x^n+m^  if  n  =  ’L^AfvE',  M  = 
[ :IS  _°Im  ] ,  and  ’L,  T-1  £  x  ("+m) . 

Reference  2  in  the  publications  funded  by  this  project  (Section  7)  provides  sufficient  conditions  for  the 
existence  of  a  J-spectral  factorization.  These  conditions  are  used  in  the  main  result  stated  below. 


Theorem  4  Let  n  =  n~  £  and  partition  as 

ML“.  IfUn(juj)  >  0  and  n 22{ju)  <  0  Mu  £lU{oc},  then 


where  Hu  £  and  YI22  £ 


1.  n  has  a  Jn,m- spectral  factorization  (T,  M). 

2.  The  Jn^m- spectral  factorization  (T,Af)  is  a  hard  factorization  of  n. 

3.  For  G  £  let  (A,B,C,D)  denote  the  state-space  realization  o/4'[<f]  in  Equation  5.  All 

solutions  P  =  P1  to  the  KYP  LMI  (Equation  6)  satisfy  P  >  0. 


Factorization  conditions  in  [4,  5]  connect  classical  passivity  multipliers  and  their  IQC  counterparts. 
Theorem  4  provides  a  connection  between  classical  passivity  multipliers  and  dissipation  theory.  Specifically, 
let  H  be  a  classical  passivity  multiplier  proving  stability  for  the  interconnection  of  G  and  a  finite-gain  system 

A.  It  follows  by  a  simple  perturbation  argument,  e.g.  as  in  [4],  that  stability  can  be  demonstrated  with  the 

eI  U*  1 

fr  -e  j  .  The  conditions  in  Theorem  4  hold  for  this  multiplier 

JlAp  J 

and  thus  a  J-spectral  factorization  of  n  exists.  Moreover,  there  is  a  dissipation  inequality  that  proves 
stability  of  the  feedback  interconnection.  In  other  words,  if  stability  can  be  demonstrated  by  a  classical 
passivity  multiplier  then  it  can  also  be  demonstrated  via  a  dissipation  inequality. 


(frequency-domain)  IQC  test  using  n  = 


3  Robustness  Analysis  for  Linear  Parameter  Varying  Systems 


The  main  result  in  the  previous  section  connects  dissipation  inequalities  and  integral  quadratic  constraints. 
As  mentioned  previously,  this  enables  new  applications  of  IQCs  to  analyze  the  robustness  of  time-varying 
and  nonlinear  systems.  This  section  considers  the  robustness  of  uncertain  linear  parameter  varying  (LPV) 
systems.  Details  on  the  results  contained  in  this  section  can  be  found  in  Reference  3  of  the  publications 
generated  by  this  research  (Section  7). 

The  uncertain  system  is  described  by  the  feedback  interconnection  of  an  LPV  system  G  and  an  uncer¬ 
tainty  A.  This  feedback  interconnection  with  A  wrapped  around  the  top  of  G  is  denoted  FU(G,  A).  The 
LPV  system  G  is  a  linear  system  whose  state  space  matrices  depend  on  a  time-varying  parameter  vector 
p  :  M+  — >  as  follows: 


xG(t)  =  AG{p(t ))  xG(t)  +  BG(p{t )) 


w{t.) 

d{t) 


v(t) 

e(t) 


CG(p(t))  xG(t)  +  DG(p(t)) 


■ w{t ) 
d{t ) 


(21) 


where  xG  £  MnG  is  the  state,  w  £  and  d  £  Mnd  are  inputs,  and  v  £  Wlv  and  e  £  Mn®  are  outputs.  The 
state  matrices  of  G  have  dimensions  compatible  with  these  signals,  e.g.  AG(p)  £  MnGxnG_  jn  addition,  the 
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state  matrices  are  assumed  to  be  continuous  functions  of  p.  The  state  matrices  at  time  t  depend  on  the 
parameter  vector  at  time  t.  Hence,  LPV  systems  represent  a  special  class  of  time-varying  systems.  The 
explicit  dependence  on  t  is  occasionally  suppressed  to  shorten  the  notation.  Moreover,  it  is  important  to 
emphasize  that  the  state  matrices  are  allowed  to  have  an  arbitrary  dependence  on  the  parameters.  This 
is  called  a  “gridded”  LPV  system  and  is  more  general  than  “LFT”  LPV  systems  whose  state  matrices  are 
restricted  to  have  a  rational  dependence  on  the  parameters  [1, 12, 15]. 

The  parameter  p  is  assumed  to  be  a  continuously  differentiable  function  of  time  and  admissible  tra¬ 
jectories  are  restricted  to  a  known  compact  set  V  C  Wlp.  In  addition,  the  parameter  rates  of  variation 
p  :  M+  — »•  V  are  assumed  to  lie  within  a  hyperrectangle  V  :=  {q  E  <  qt  <  Lj,  i  =  1, . . .  ,np}.  The 

set  of  admissible  trajectories  is  defined  as 

T  :=  j/9  :  M+  ->  Rnp  :  p  G  C1,  p(t)  €  V  and  p{t)  E  V  Vt  >  oj  (22) 

The  parameter  trajectory  is  said  to  be  rate  unbounded  if  V  =  Mnp. 

Throughout  the  section  it  is  assumed  that  the  uncertain  system  has  a  form  of  nominal  stability.  Specif¬ 
ically,  G  is  assumed  to  be  parametrically-dependent  stable  as  defined  in  [25]. 

Definition  5  G  is  parametrically-dependent  stable  if  there  is  a  continuously  differentiable  function  P  : 
Mnp  — >■  snGXnc  such  that  P(p)  >  0  and 

np  BP 

AG(P)TP(P)  +  P{p)AG(p)  +  g-Qi  <  0  (23) 

i=i  °Pl 

hold  for  all  p  G  P  and  all  q  G  V . 

As  discussed  in  [25],  parametric-stability  implies  G  has  a  strong  form  of  robustness.  In  particular,  the 
state  xcif)  of  the  autonomous  response  (w  =  0,  d  =  0)  decays  exponentially  to  zero  for  any  initial 
condition  xq{ 0)  G  Mn<3  and  allowable  trajectory  p  G  T  (Lemma  3.2.2  of  [25]).  Moreover,  the  state  xc{t) 
of  the  forced  response  decays  asympotically  to  zero  for  any  initial  condition  xc(0)  G  Mn<3,  allowable 
trajectory  p  E  T,  and  inputs  w,d  G  L-2  (Lemma  3.3.2  of  [24]).  The  parameter-dependent  Lyapunov 
function  V(xg,p)  ■=  XqP(p)xg  plays  a  key  role  in  the  proof  of  these  results.  To  shorten  the  notation,  a 
differential  operator  BP  :  V  x  V  — >  Mna:  is  introduced  as  in  [17].  BP  is  defined  as  8P(p ,  q)  :=  X^=i  §jr(p)  (h ■ 
This  simplifies  the  expression  of  Lyapunov-type  inequalities  similar  to  Equation  23. 

The  uncertainty  A  :  Lff  [0,  oo)  — >  Vff  [0,  oo)  is  a  bounded,  causal  operator.  The  notation  A  is  used 
to  denote  the  set  of  bounded,  causal  uncertainties  A.  The  input/output  behavior  of  the  uncertain  set  is 
bounded  using  quadratic  constraints  as  described  further  in  the  next  section.  At  this  point  it  is  sufficient 
to  state  that  A  can  have  block-structure  as  is  standard  in  robust  control  modeling  [27].  A  can  include 
blocks  that  are  hard  nonlinearities  (e.g.  saturations)  and  infinite  dimensional  operators  (e.g.  time  delays) 
in  addition  to  true  system  uncertainties.  The  term  uncertainty  is  used  for  simplicity  when  referring  to  the 
perturbation  A. 

The  objective  of  this  section  is  to  assess  the  robustness  of  the  uncertain  system  FU(G ,  A).  For  a  given 
A  G  A,  the  induced  L2  gain  from  d  to  e  is  defined  as: 

\\Pu(G,  A)  ||  :=  sup  (24) 

0^deL2d[0,oo)  H'hl2 
peT,  £g(0)=0 

Two  forms  of  robustness  are  considered.  First,  the  worst-case  induced  L2  gain  from  input  d  to  the  output 
e  is  defined  as 


sup  \\FU(GP,  A)  || .  (25) 

AeA 
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This  is  the  worst-case  gain  over  all  uncertainties  A  G  A  and  admissible  trajectories  p  G  T.  Second, 
the  system  has  robust  asymptotic  stability  if  xc(t)  — *  0  for  any  initial  condition  xq{ 0)  G  Mn<3,  allowable 
trajectory  p  G  T,  disturbance  d  G  L-2  and  uncertainty  A  G  A. 

The  main  result  (Theorem  5  below)  provides  a  sufficient  condition  for  when  an  uncertain  LPV  system 
has  both  robust  asymptotic  stability  and  bounded  worst-case  gain.  The  results  is  based  on  the  intercon¬ 
nection  of  the  nominal  LPV  system  G  and  the  filter  T  of  the  IQC  factorization,  similar  to  Figure  2.  The 
dynamics  of  this  interconnection  are  described  by  w  =  A(w)  and 

x  =  A(p)x  +  Bi(p)w  +  B2(p)d 

z  =  Ci(p)x  +  Du(p)w  +  D12(p)d  (26) 

e  =  C2{p)x  +  D2i(p)w  +  D22(p)d, 

with  x  :=  \Xp]  G  MnG+nv>  is  the  extended  state.  Removing  the  uncertainty  A  from  the  analysis  intercon¬ 
nection,  w  can  be  viewed  as  an  external  signal  subject  to  the  constraint  /Qr  zT (t)M z(t)  dt  >  0. 


Theorem  5  Let  G  be  a  parametrically  stable  LPV  system  defined  by  eq.  (21)  and  A  :  Ltf"  [0,  oc)  — x 
0,  °o)  be  a  bounded,  causal  operator  such  that  FU(G,  A)  is  well-posed.  Assume  A  satisfies  the  IQC 
parameterized  by  11(A)  =  \P~M( A)T  with  'P  stable.  If 


1. 


The  combined  multiplier,  partitioned  as  n(A)  = 
Vcj  G  M  U  {oo}  where  Iln  is  nv  x  nv  and  II22  is  n 


riii  1I12 

n~2  1I22  ’ 

w  x  nw . 


satisfies  II 1 1  (joj)  >  0  and  TL22  (juj)  <  0 


2.  There  exists  a  continuously  differentiable  P  :  V  — >  Sn'xX1lx ,  and  a  scalar  7  >  0  such  that 


’  ATP+PA+dP  PBi  PB2  " 

cf 

\cf-\ 

T 

r  ct  1 

~cf- 

BjP  0  0 

BfP  0  -72/ 

+ 

DL 

DL_ 

to  to 

1 _ 

+ 

Dfi 

M  (A) 

Ofi 

DL. 

hold  for  all  p  G  V  and  all  q  G? . 

Then 

a)  For  any  x(0)  G  MnG+nv>  and  d  G  L2,  limj'_5.00  x(T)  =  0 


(27) 


b)  \\FU(G,  A)  ||  <7 


The  implementation  of  Theorem  5  involves  some  numerical  issues.  These  are  briefly  described  here. 
If  the  IQC  is  parameterized  such  that  M( A)  is  an  affine  function  of  A  then  theorem  5  involves  parameter 
dependent  LMI  conditions  in  the  variables  P(p)  and  A.  A  needs  to  satisfy  condition  1  in  theorem  5,  i.e. 
II11  >  0  and  II22  <  0.  These  are  infinite  dimensional  (one  LMI  for  each  (p,  q )  GfxP)  and  they  are 
typically  approximated  with  finite-dimensional  LMIs  evaluated  on  a  grid  of  parameter  values.  Additionally, 
the  main  decision  variable  is  the  function  P{p)  which  must  be  restricted  to  a  finite  dimensional  subspace. 
A  common  practice  [2,26]  is  to  restrict  P(p)  to  be  a  linear  combination  of  user-specified  basis  functions. 
The  analysis  can  then  be  performed  as  a  finite-dimensional  SDP  [3],  e.g.  minimizing  7  subject  to  the 
approximate  finite-dimensional  LMI  conditions.  This  paper  focused  on  gridded  LPV  systems  whose  state 
matrices  have  an  arbitrary  dependence  on  the  parameter.  If  the  LPV  system  has  a  rational  dependence  on 
p  then  finite  dimensional  LMI  conditions  can  be  derived  (with  no  gridding)  using  the  techniques  in  [1, 12]. 
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4  Nonlinear  Robustness  Analysis 


The  main  result  in  Section  2  connects  dissipation  inequalities  and  integral  quadratic  constraints.  As 
mentioned  previously,  this  enables  new  applications  of  IQCs  to  analyze  the  robustness  of  time- varying  and 
nonlinear  systems.  This  section  considers  the  analysis  of  nonlinear  systems.  Details  on  the  results  contained 
in  this  section  can  be  found  in  Reference  1  of  the  publications  generated  by  this  research  (Section  7). 

Consider  a  nonlinear  system  governed  by  differential  equations  of  the  follwoing  form 

x{t)  =  f(x(t ))  +  g(x(t))w(t), 
z(t)  =  h(x(t)), 

where  ieR,  x(0)  =  xq  G  Mn,  x(t)  G  Mn,  z(t)  G  Rp,  w(t)  G  Mm.  The  functions  /  :  Mn  — >  Mn,g  :  Mn  — >  Mnxm 
and  h  :  Mn  — >  MP  are  assumed  to  be  Lipschitz  continuous,  or  locally  Lipschitz  continuous,  depending  on 
the  situation.  If  /  and  g  are  not  Lipschitz  continuous  (as  in  the  case  of  polynomial  /  and  g,  for  example), 
then  the  differential  equation  may  exhibit  finite  escape  times  in  the  presence  of  bounded  inputs  and/or 
initial  conditions. 

The  goal  of  this  research  is  quantitative,  local  analysis  of  nonlinear  dynamical  systems.  By  “quantita¬ 
tive”  we  mean  algorithms  and  sufficient  conditions  which  lead  to  concrete  guarantees  about  a  particular 
system’s  response.  By  “local”  we  refer  to  guarantees  about  the  reachability  and/or  system  gain  which  are 
predicated  on  assumptions  concerning  the  magnitude  of  initial  conditions  and  input  signals.  We  extensively 
use  the  basic,  fundamental  ideas  from  dissipative  systems  theory  [22],  [7],  barrier  functions  and  reachabil¬ 
ity  [19],  [13],  and  nonlinear  optimal  control  [18],  [6].  Specifically,  we  employ  inequalities  involving  the  Lie 
derivative  of  a  scalar  function,  the  storage  function,  that  hold  throughout  regions  of  the  state  and  input 
space,  which  when  integrated  over  trajectories  of  the  system,  give  certificates  of  input/output  properties 
of  the  system.  The  necessity  of  the  existence  of  such  storage  functions  to  prove  input/output  properties, 
which  leads  to  the  most  elegant  results  of  the  above  mentioned  works,  is  actually  not  used  in  this  paper. 
Our  computational  approach  is  based  on  polynomial  storage  functions  of  fixed  degree  which  can  be  viewed 
as  extensions  of  known  linear  matrix  inequality  conditions  to  compute  reachable  sets  and  input/output 
gains  for  linear  systems  [3]. 

The  contributions  are  as  follows:  a  dissipation  inequality  formulation  of  local  reachability  and  dissi¬ 
pativeness  for  uncertain  systems  that  are  not  nominally  globally  stable  are  derived;  refinements  on  the 
reachability  and  L2  gain  conditions  that  can  be  used  to  efficiently  compute  improved  quantitative  perfor¬ 
mance  bounds;  sum-of-squares  (SOS)  characterizations  of  the  required  set  containment  conditions  in  the 
dissipation  inequalities;  proof  of  guaranteed  feasability  of  the  SOS  conditions  for  systems  with  stable  lin¬ 
earizations;  development  of  a  scheme  to  find  feasible  solutions  to  the  bilinear  SOS  conditions,  and  improve 
the  objective  through  a  specific  iteration  scheme;  and  a  collection  of  illustrative  and  realistic  examples 
illustrating  the  methods. 

One  result  on  reachability  is  provided  to  demonstrate  the  basic  approach.  Details  on  the  remaining 
results  can  be  found  in  Reference  1  of  the  publications  generated  by  this  research  (Section  7).  Specifically, 
we  establish  conditions  which  guarantee  invariance  of  certain  sets  under  L2  and  pointwise-in-time  (Loo- 
like)  constraints  on  w.  These  are  subsequently  referred  to  as  “reachability”  results,  since  the  conclusions 
yield  outer  bounds  on  the  set  of  reachable  states.  In  that  vein,  w  is  interpreted  as  a  disturbance,  whose 
worst-case  effect  on  the  state  x  is  being  quantified.  We  obtain  bounds  on  x  that  are  tightly  linked  with 
the  assumed  bounds  on  w  and  xo,  and  specifically  allow  for  systems  which  are  not  well-defined  on  all  input 
signals  (finite  escape  times). 

A  known  set  VV  C  Mm  is  used  to  express  any  Loo-like,  pointwise-in-time  bound  on  the  signal  w,  namely 
w(t )  G  W  for  all  t.  Setting  VV  =  Rm  is  equivalent  to  the  absence  of  known,  pointwise-in-time  bounds  on 
w. 


DISTRIBUTION  A:  Distribution  approved  for  public  release. 


Theorem  6  Suppose  W  C  Mm.  Assume  that  f  and  g  in  (28)  are  Lipschitz  continuous  on  Mn.  Suppose 
t  >  0,  and  a  differentiable  Q  :  Mn  — >  M  satisfies  Q( 0)  <  r2  and 

2  x  W  C  |(.x,  w)  G  Mn  x  Mm  :  VQ(x)  ■  [f(x)  +  g(x)w]  <  wTw }  .  (29) 

Consider  xo  G  W^1  Q(x o)  <  t2  and  w  G  L™  with  w(t)  G  W  for  all  t.  If  ||u;||2  <  r2  —  Q(x o),  the 

solution  to  (28)  with  x(0)  =  xq  satisfies  Q(x(t))  <  r2  for  all  t,  and  hence  x(t)  G  Qq’%  for  all  t. 

Without  loss  of  generality,  Q  in  Theorem  6  can  be  taken  to  be  zero  at  x  =  (1.  For  instance,  define 
Q(x)  :=  Q(x)  —  Q{ 0)  and  f2  :=  r2  —  Q(0).  The  conditions  of  Theorem  6  hold  with  Q  replacing  Q,  and  the 
same  norm  bound  (i.e.  reachable  set)  is  obtained.  Computational  approaches  based  on  the  S-procedure 
and  sum-of-squares  are  described  further  in  the  paper. 

5  IQC  Analysis  for  Certifying  Exponential  Convergence 

The  standard  stability  result  for  interconnections  involving  nonlinearities  described  by  IQCs  [9-11]  provides 
a  frequency-domain  test  for  certifying  BIBO  stability  of  the  associated  interconnected  system.  Since  the 
frequency-domain  test  is  difficult  to  verify  computationally,  a  standard  approach  is  to  use  the  Kalman- 
Yakubovich-Popov  (KYP)  lemma  to  obtain  an  equivalent  linear  matrix  inequality  (LMI)  that  can  then  be 
checked  using  a  conventional  convex  programming  software  package. 

The  approach  outlined  above  is  useful  in  verifying  the  robust  stability  of  an  interconnected  system 
containing  components  that  are  nonlinear,  uncertain,  poorly  modeled,  or  otherwise  problematic.  With 
minor  modifications,  the  same  result  can  be  adapted  for  use  in  certifying  L2  gain  bounds  or  passivity. 
Unfortunately,  the  approach  cannot  prove  exponential  stability.  To  contrast  both  notions  of  stability,  we 
have  (in  discrete  time): 


L2  stability  implies: 
exponential  stability  implies: 


OO 

||.Tfc||2  <  y2  for  some  7,  while 

k= 0 

||xfc||  <  cXk  for  some  c  >  0  and  A  G  (0, 1) 


Note  that  in  the  L2  case,  the  norm  is  invariant  under  rearrangement.  Thus,  two  signals  with  the  same  7 
both  eventually  go  to  zero,  but  their  transient  behaviors  may  be  completely  different.  In  the  exponential 
case,  however,  two  signals  with  the  same  (c,  A)  have  the  same  decaying  exponential  envelope.  Therefore, 
exponential  stability  is  much  stronger  than  its  L2  counterpart,  and  it  is  useful  to  be  able  to  certify  it. 

The  standard  IQC  result  can  only  certify  L2  stability.  The  authors  of  [Megretski  and  Rantzer  ’97] 
address  this  shortcoming  by  pointing  out  that  BIBO  stability  often  implies  exponential  stability  in  many 
cases  of  practical  interest.  So  by  certifying  an  L2  gain  bound,  we  often  get  exponential  stability  for  free. 
However,  in  such  cases,  one  can  only  prove  the  existence  of  some  (c,  A).  One  could  not,  for  example, 
optimize  to  find  the  certificate  with  the  fastest  possible  exponential  rate  (i.e.  minimizing  A). 

Under  the  support  of  AFOSR,  we  developed  an  IQC-based  method  for  certifying  exponential  conver¬ 
gence  that  also  allows  one  to  optimize  over  A.  We  also  provide  nontrivial  instances  for  which  our  exponential 
rate  bounds  are  tight. 

To  illustrate  our  approach,  it  is  useful  to  first  consider  a  simple  case.  Consider  a  discrete  linear  time- 
invariant  (LTI)  plant  G  with  state-space  realization  (A,  B,  C ,  D ).  Suppose  G  is  connected  in  feedback  with 
a  passive  nonlinearity  A.  A  sufficient  condition  for  BIBO  stability  is  that  there  exists  a  positive  definite 
matrix  P  y  0  and  a  scalar  A  >  0  satisfying  the  linear  matrix  inequality  (LMI) 
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If  we  define  V(x)  :=  xTPx,  then  (30)  implies  that  V  decreases  along  trajectories:  V(xk+{)  <  V{xk )  for 
all  k.  BIBO  stability  then  follows  from  positivity  and  boundedness  of  V.  But  observe  that  when  (30) 
holds,  we  may  replace  the  right-hand  side  by  —sP  for  some  e  >  0  sufficiently  small.  We  then  conclude 
that  V{xk+\)  <  (1  —  e)V(xk)  for  all  k  and  exponential  stability  follows.  We  can  then  maximize  e  subject 
to  feasibility  of  (30)  to  further  improve  the  rate  bound.  A  simple  way  of  carrying  out  this  maximization 
is  to  perform  a  bisection  search  on  e,  since  (30)  is  a  convex  program  for  every  fixed  e. 

The  simple  trick  shown  above  works  because  passivity  is  a  very  special  case  of  an  IQC  for  which  there 
are  no  dynamics  involved.  In  other  words,  passivity  can  be  verified  in  a  pointwise  fashion  by  checking  that 
uTyk  >  0  for  all  k.  Unfortunately,  the  trick  shown  above  does  not  work  in  the  general  dynamic  IQC  setting 
due  to  the  different  role  played  by  P  in  the  associated  LMI.  The  LMI  used  in  IQC  theory  comes  from  the 
KYP  lemma  and  although  it  is  structurally  similar  to  (30),  P  is  not  positive  definite  in  general  and  V  may 
not  decrease  along  trajectories. 

Our  key  insight  is  that  with  a  suitable  modification  to  both  the  LMI  and  the  IQC  definition,  we  obtain 
a  condition  that  can  certify  exponential  stability  for  a  given  rate  A.  The  resulting  optimization  problem  is 
convex  for  any  A,  so  we  can  optimize  over  A  by  performing  a  bisection  search.  We  show  that  our  modified 
IQC  approach  finds  tight  exponential  bounds  for  a  simple  yet  nontrivial  example:  a  third-order  plant  in 
feedback  with  a  saturating  nonlinearity.  Additional  details  on  the  results  contained  in  this  section  can  be 
found  in  Reference  4  of  the  publications  generated  by  this  research  (Section  7). 


6  Analyzing  Optimization  Algorithms 


Electromechanical  systems  often  contain  embedded  optimization  algorithms.  One  example  is  modern 
applications  of  model-predictive  control.  Other  examples  include  machine  learning  and  computer  vision, 
where  large  quantities  of  data  must  be  processed  in  real-time  on  a  small  platform  such  as  a  robot  or  UAV. 
In  all  such  cases,  the  optimization  algorithms  are  sequential  in  nature,  and  are  carried  out  either  for  a  fixed 
number  of  steps,  or  until  a  desired  error  tolerance  is  reached. 

Sequential  algorithms  can  be  interpreted  as  uncertain  dynamical  systems,  where  the  uncertainty  is  due 
to  the  function  being  optimized,  and  the  dynamics  are  due  to  the  choice  of  sequential  algorithm.  As  a  simple 
example,  consider  the  heavy  ball  method  applied  to  some  unknown  differentiable  function  /  :  Mn  — *•  M. 

xk+x  =  xk  -  aS7f(xk)  +  /3(xk  -  xfc_i)  (31) 


Here,  a  and  f3  are  the  stepsize  and  momentum  parameters,  respectively.  By  including  additional  definitions 
for  pk ,  ukl  yk,  we  may  rewrite  (31)  as  follows: 
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If  we  let  G  be  the  linear  time-invariant  system  G  :  u  y  described  by  (32a)-(32b)  and  we  let  V/  be  the 
static  map  y  u  described  by  (32c),  then  we  can  represent  (31)  as  the  block  diagram  of  Figure  3. 

It  turns  out  that  many  different  kinds  of  iterative  algorithms  may  be  abstracted  in  the  form  of  Figure  3. 
These  include  gradient-based  schemes  such  as  gradient  descent  and  its  accelerated  variants  such  as  the 
heavy  ball  method  and  the  fast  gradient  method  (also  known  as  Nesterov’s  accelerated  method).  More 
complicated  examples  include  projected  variants  for  use  in  constrained  optimization,  proximal  algorithms, 
and  operator-splitting  methods  such  as  the  alternating  direction  method  of  multipliers  (ADMM). 

A  common  way  of  evaluating  the  performance  of  an  iterative  optimization  algorithm  is  to  bound  its 
worst-case  convergence  rate.  When  viewed  as  a  dynamical  system,  this  amounts  to  proving  exponential 
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Figure  3:  Block  diagram  representing  the  uncertain  dynamics  of  the  heavy  ball  method  (31).  G  represents 
the  LTI  dynamics  of  the  iterative  algorithm  while  V/  is  the  gradient  of  the  function  being  optimized. 

stability  and  computing  the  smallest  associated  decay  rate  A  G  (0, 1).  Therefore,  the  method  described  in 
Section  5  is  perfectly  suited  for  this  analysis. 

When  we  carry  out  the  analysis,  we  recover  the  best  possible  rates  for  a  multitude  of  algorithms  as 
found  in  the  existing  literature  on  optimization  algorithms.  Moreover,  the  same  proof  works  for  all  the 
algorithms.  In  a  further  study,  we  used  the  same  IQC-based  approach  to  analyze  the  ADMM  algorithm 
specifically.  We  found  provably  tight  bounds  on  its  performance  as  well  as  a  principled  approach  to 
hyperparameter  selection.  The  benefits  of  our  analysis  are  the  inherent  simplicity  of  the  approach  and  the 
fact  that  the  associated  LMIs  are  very  small  and  can  be  solved  in  milliseconds  using  generic  solvers  on 
cheap  hardware.  Additional  details  on  the  results  contained  in  this  section  can  be  found  in  References  5 
and  6  of  the  publications  generated  by  this  research  (Section  7). 

7  Publications 

The  following  is  a  list  of  publications  that  have  either  been  partially  supported  or  fully  supported  during 
the  this  AFOSR  project. 

1.  E.  Summers,  A.  Chakraborty,  W.  Tan,  U.  Topcu,  P.  Seiler,  G.  Balas,  and  A.  Packard,  “Quantitative 
local  L2-gain  and  Reachability  analysis  for  nonlinear  systems,”  Int.  Journal  of  Robust  and  Nonlinear 
Control,  vol.  23,  p.  1115-1135,  2013. 

2.  P.  Seiler,  “Stability  Analysis  with  Dissipation  Inequalities  and  Integral  Quadratic  Constraints,”  to 
appear  in  the  IEEE  Transactions  on  Automatic  Control,  2015. 

3.  H.  Pfifer  and  P.  Seiler,  “Less  Conservative  Robustness  Analysis  of  Linear  Parameter  Varying  Sys¬ 
tems  Using  Integral  Quadratic  Constraints,”  submitted  to  the  International  Journal  of  Robust  and 
Nonlinear  Control,  2015. 

4.  R.  Boczar,  L.  Lessard,  and  B.  Recht,  “Exponential  convergence  bounds  using  integral  quadratic 
constraints,”  submitted  to  the  2015  IEEE  Conference  on  Decision  and  Control. 

5.  R.  Nishihara,  L.  Lessard,  B.  Recht,  A.  Packard,  and  M.  I.  Jordan,  “A  general  analysis  of  the  conver¬ 
gence  of  ADMM,”  to  appear  in  the  2015  International  Conference  on  Machine  Learning. 

6.  L.  Lessard,  B.  Recht,  and  A.  Packard,  “Analysis  and  design  of  optimization  algorithms  via  integral 
quadratic  constraints,”  submitted  to  the  SIAM  Journal  on  Optimization,  2015. 

References 

[1]  P.  Apkarian  and  P.  Gahinet.  A  convex  characterization  of  gain-scheduled  H0 0  controllers.  IEEE  Trans, 
on  Automatic  Control,  40:853-864,  1995. 


DISTRIBUTION  A:  Distribution  approved  for  public  release. 


[2]  G.  Balas.  Linear,  parameter-varying  control  and  its  application  to  a  turbofan  engine.  Int.  Journal  of 
Robust  and  Nonlinear  Control ,  12:763-796,  2002. 

[3]  S.  Boyd,  L.  El  Ghaoui,  E.  Feron,  and  V.  Balakrishnan.  Linear  Matrix  Inequalities  in  System  and 
Control  Theory ,  volume  15  of  Studies  in  Applied  Mathematics.  SIAM,  1994. 

[4]  J.  Carrasco,  W.P.  Heath,  and  A.  Lanzon.  Factorization  of  multipliers  in  passivity  and  IQC  analysis. 
Automatica,  48(5):909-916,  2012. 

[5]  M.  Fu,  S.  Dasgupta,  and  Y.C.  Soh.  Integral  quadratic  constraint  approach  vs.  multiplier  approach. 
Automatica ,  41:281-287,  2005. 

[6]  J.W.  Helton  and  M.  James.  Extending  H0 0  control  to  nonlinear  systems:  control  of  nonlinear  systems 
to  achieve  performance  objectives.  Frontiers  in  Applied  Mathematics.  SIAM,  1999. 

[7]  D.J.  Hill  and  P.J.  Moylan.  Dissipative  dynamical  systems:  Basic  input-output  and  state  properties. 
Journal  of  the  Franklin  Institute ,  309(5)  :327-357,  1980. 

[8]  H.K.  Khalil.  Nonlinear  Systems.  Prentice  Hall,  third  edition,  2001. 

[9]  A.  Megretski  and  A.  Rantzer.  System  analysis  via  integral  quadratic  constraints:  Part  I.  Technical 
Report  TFRT-7531-SE,  Lund  Inst,  of  Technology,  1995. 

[10]  A.  Megretski  and  A.  Rantzer.  System  analysis  via  integral  quadratic  constraints:  Part  II.  Technical 
Report  TFRT-7559-SE,  Lund  Inst,  of  Technology,  1995. 

[11]  A.  Megretski  and  A.  Rantzer.  System  analysis  via  integral  quadratic  constraints.  IEEE  Trans,  on 
Aut.  Control,  42(6):819-830,  1997. 

[12]  A.  Packard.  Gain  scheduling  via  linear  fractional  transformations.  Systems  and  Control  Letters , 
22:79-92,  1994. 

[13]  S.  Prajna.  Barrier  certificates  for  nonlinear  model  validation.  Automatica ,  42(1):117-126,  2006. 

[14]  A.  Rantzer.  On  the  Kalman-Yakubovich-Popov  lemma.  Systems  and  Control  Letters ,  28(1):7— 10, 
1996. 

[15]  C.  Scherer.  Advances  in  linear  matrix  inequality  methods  in  control ,  chapter  Robust  mixed  control 
and  linear  parameter-varying  control  with  full-block  scalings,  pages  187-207.  SIAM,  2000. 

[16]  C.  Scherer  and  S.  Weiland.  Linear  matrix  inequalities  in  control,  2000.  Version  3.0. 

[17]  C.  Scherer  and  S.  Wieland.  Linear  matrix  inequalities  in  control.  Lecture  notes  for  a  course  of  the 
Dutch  institute  of  systems  and  control,  Delft  University  of  Technology,  2004. 

[18]  A.  J.  van  der  Schaft.  C2-gain  and  Passivity  Techniques  in  Nonlinear  Control.  Springer- Verlag,  New 
York  and  Berlin,  second  edition,  2000. 

[19]  R.  Vinter.  A  characterization  of  the  reachable  set  for  nonlinear  control  systems.  SIAM  Journal  on 
Control  and  Optimization ,  18(6):599-610,  1980. 

[20]  J.C.  Willems.  Least  squares  stationary  optimal  control  and  the  algebraic  Riccati  equation.  IEEE 
Trans,  on  Aut.  Control ,  16:621-634,  1971. 

[21]  J.C  Willems.  Dissipative  dynamical  systems  part  i:  General  theory.  Archive  for  Rational  Mech.  and 
Analysis ,  45(5):321-351,  1972. 


DISTRIBUTION  A:  Distribution  approved  for  public  release. 


[22]  J.C.  Willems.  Dissipative  dynamical  systems  part  I:  General  theory.  Archive  for  Rational  Mechanics 
and  Analysis ,  45(5):321-351,  1972. 

[23]  J.C  Willems.  Dissipative  dynamical  systems  part  ii:  Linear  systems  with  quadratic  supply  rates. 
Archive  for  Rational  Mech.  and  Analysis,  45(5):352-393,  1972. 

[24]  F.  Wu.  Control  of  Linear  Parameter  Varying  Systems.  PhD  thesis,  University  of  California,  Berkeley, 
1995. 

[25]  F.  Wu.  Control  of  Linear  Parameter  Varying  Systems.  PhD  thesis,  University  of  California,  Berkeley, 
1995. 

[26]  F.  Wu,  X.  H.  Yang,  A.  Packard,  and  G.  Becker.  Induced  C-2  norm  control  for  LPV  systems  with 
bounded  parameter  variation  rates.  Int.  Journal  of  Robust  and  Nonlinear  Control,  6:983-998,  1996. 

[27]  K.  Zhou,  J.C.  Doyle,  and  K.  Glover.  Robust  and  Optimal  Control.  Prentice-Hall,  1996. 


DISTRIBUTION  A:  Distribution  approved  for  public  release. 


Response  ID:4606  Data 


1. 

1.  Report  Type 
Final  Report 
Primary  Contact  E-mail 

Contact  email  if  there  is  a  problem  with  the  report. 

seileOI  7@umn.edu 

Primary  Contact  Phone  Number 

Contact  phone  number  if  there  is  a  problem  with  the  report 

612-626-5289 

Organization  /  Institution  name 

University  of  Minnesota 

Grant/Contract  Title 

The  full  title  of  the  funded  effort. 

A  Merged  IQC/SOS  Theory  for  Analysis  and  Synthesis  of  Nonlinear  Control  System 

Grant/Contract  Number 

AFOSR  assigned  control  number.  It  must  begin  with  "FA9550"  or  "F49620"  or  "FA2386". 

FA9550-1 2-1 -0339 

Principal  Investigator  Name 

The  full  name  of  the  principal  investigator  on  the  grant  or  contract. 

Peter  Seiler 

Program  Manager 

The  AFOSR  Program  Manager  currently  assigned  to  the  award 

Fariba  Fahroo 

Reporting  Period  Start  Date 

06/1 5/2012 

Reporting  Period  End  Date 

06/1 4/2015 

Abstract 

The  connections  between  dissipation  inequalities  and  integral  quadratic  constraints  (IQCs)  was  the  major 
thrust  explored  during  the  program.  In  particular,  it  was  shown  that  existing  frequency  domain  IQC  analysis 
conditions  are  equivalent  (under  mild  technical  assumptions)  to  a  related  time-domain  dissipation 
inequality  condition.  This  time-domain  approach  enables  applications  of  IQCs  to  analyze  the  robustness  of 
uncertain  time-varying  and  nonlinear  systems.  These  theoretical  techniques  can  be  used  to  improve  the 
design  and  robustness  of  advanced  flight  control  algorithms. 

Distribution  Statement 

This  is  block  12  on  the  SF298  form. 

Distribution  A  -  Approved  for  Public  Release 

Explanation  for  Distribution  Statement 

If  this  is  not  approved  for  public  release,  please  provide  a  short  explanation.  E.g.,  contains  proprietary  information. 
SF298  Form 

Please  attach  your  SF298  form.  A  blank  SF298  can  be  found  here.  Please  do  not  password  protect  or  secure  the  PDF 
The  maximum  file  size  for  an  SF298  is  50MB. 

%SS¥rIBUTION  A:  Distribution  approved  for  public  release. 


Upload  the  Report  Document.  File  must  be  a  PDF.  Please  do  not  password  protect  or  secure  the  PDF  .  The 
maximum  file  size  for  the  Report  Document  is  50MB. 

AFOSR_FinalReport.pdf 

Upload  a  Report  Document,  if  any.  The  maximum  file  size  for  the  Report  Document  is  50MB. 

Archival  Publications  (published)  during  reporting  period: 

Changes  in  research  objectives  (if  any): 

Change  in  AFOSR  Program  Manager,  if  any: 

Extensions  granted  or  milestones  slipped,  if  any: 

AFOSR  LRIR  Number 
LRIR  Title 
Reporting  Period 
Laboratory  Task  Manager 
Program  Officer 
Research  Objectives 
Technical  Summary 

Funding  Summary  by  Cost  Category  (by  FY,  $K) 


Starting  FY 

FY+1 

FY+2 

Salary 

Equipment/Facilities 

Supplies 

Total 

Report  Document 

Report  Document  -  Text  Analysis 

Report  Document  -  Text  Analysis 

Appendix  Documents 

2.  Thank  You 

E-mail  user 

May  28,  201 5  12:32:04  Success:  Email  Sent  to:  seileOI  7@umn.edu 


DISTRIBUTION  A:  Distribution  approved  for  public  release. 


